Calculation of Magnetic Penetration Depth Length A(T) in High Tc Superconductors 
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The notion of a finite pairing interaction energy range via Nam, results in the incomplete conden- 
sation in which not all states participate in pairings. The states not participating in pairings are 
shown to yield the low energy states responsible for the linear T dependence of superelectron density 
at low T in a s-wave superconductor. We present extensive quantitative calculations of A(T) for 
all T ranges, in good agreements with experiments. It is not necessary to have nodes in the order 
parameter, to account for the linear T dependence of A(T) at low T in high Tc superconductors. 
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One of crucial parameters in a superconductor is the 
magnetic penetration depth length A(T) which reflects 
the condensation carrier density, superelectron density 
p s (T), in the London model as 



Ps(T)/pM = [X(0)/X(T)Y 



(1) 



The p s (T) plays an important role for understanding 
the nature of condensation. In the Gorter and Casimir 
two fluid model (GC), p s (T) varies as 1 - (T/T c ) 4 . But 
the BCS-p s (T) has an activation form at low T via the 
order parameter A which indicates the excitation energy 
gap. The measurements Q-o of A(T) at low T in high 
T c superconductors (HTS) are compatible with neither 
the BCS result nor the GC picture. Data indicate the 
linear T dependence of p s (T) at low T. This linear T de- 
pendence of A(T) in fact is taken as providing evidence 
that the order parameter has nodes, suggesting the d- 
wave pairing state 0. On the other hand, one of us [^]|| 
has shown that the notion of a finite pairing interaction 
energy range Td results in the incomplete condensation 
and the low energy states responsible for the linear T de- 
pendence of p s {T) at low T in a s-wave superconductor 
fL0| . Moreover, the incomplete condensation yields the 
multi-connected superconductors(MS) which can ac- 
count for the 7r-phase shift in Pb-YBCO SQUID |l2) and 
1/2 fluxoid quantum in the YBCO ring with three grain 
boundary junctions ML 

Recently, the oxygen isotope effect fl4|| , T c oc A/~ Q 
with a — 0.4 ~ 0.49 in LSCO single crystal, suggests the 
electron-phonon interaction would play an important role 
for understanding superconductivity in cuprate materi- 
als. And the BSCCO bicrystal c-asxis twist Josephson 
junction experiment |l5f| indicates the dominant order pa- 
rameter contains the s-wave and not d-wave component. 
Moreover, no node in the order parameter is observed in 
the angular dependence of the non-linear transverse mag- 
netic moment of YBCO in the Meissner state Jl6|] . On the 
other hand, the scanning tunneling microscope imaging 
the effects of individual zinc impurity atoms on super- 
conductivity in BSCCO |l7j shows the four fold symmet- 
ric quasiparticle cloud, indicating the d-wave component. 



But no four fold is observed in the same system p^] . Per- 
haps, the observation of Q may be a reflection of the 
Fermi surface. 

It is highly desirable to carry out quantitative calcula- 
tions of A(T) for all T ranges to see the accountability of 
finite Td picture for A(T) data of HTS. In this letter, we 
present extensive quantitative calculations of A(T) for all 
T ranges in good agreements with data of HTS. For this, 
it is worthy to recapitulate the pertinent results for the 
notion of a finite Td || . 

To see the phase transition, the transition tempera- 
ture T c should be a finite value, that is, neither zero nor 
infinite. To have a finite value of T c , the pairing interac- 
tion energy range Td should be finite, since T c is scaled 
with Td within the pairing theory JsJ . In other words, the 
order parameter A(fc,o>) may be written as ft)) 



A(fc,w) 



A 





for 
for 



<T d 
>T d 



(2) 



for all frequencies u>. Here ej, is the usual normal state 
excitation energy with the momentum fc, measured with 
respect to the Fermi level. 
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FIG. 1. Schematic diagram showing the order parameter 
A in the finite pairing interaction energy ranges Td view. 
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Here the natural units of Ti = c = fee = 1 are used. 
Later we use A(T) with T for A as well. Note that w has 
no constraint and that in the case of pairings of carriers 
via the electron-phonon interaction, the Td corresponds 
to the Debye temperature. However, the nature of Td in 
HTS is still unknown. Our results are not depending on 
the nature of Td- The order parameter A is the solution 
of the BCS like equation || 

1 f Td de 3 

- = / p tanh^E, 3 
9 Jo E 2 

where E — (e 2 + A 2 ) 1 / 2 ,/? = 1/T, and g corresponds to 
the BCS coupling parameter N(0)Vbcs- The solution of 
A for g are shown in the unit of Td in Fig. 0. 
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FIG. 2. The order parameter A as a function of g and T. 

The equation for T c from Eq. (|J) via A(T C ) = may 
be written as |8||jJ 

l/g=(2/n)J2W3)tan- l (y/j), (4) 

3 

where y — Td/^T C , and sum is over the positive odd 
integers j. The factor of arctangent function makes the 
sum converge. For large y, Eq. (||) yields the BCS result 
T C (BCS). The quantitative calculations of T c are given 
in Fig. |, together with the BCS T C (BCS). Unlike the 
BCS result of T C {BCS), the T c from Eq.(|) does not 
have any upper limit. The fact is that for large g > 
2.32, T c increases with increasing g as T c — gTd/2. One 
interesting value of g = 0.657 yields T c = 100 K with Td 
— 400 K which is of the order of the Debye temperature 
in HTS. This value of g may be realized in YBCO by 
considering the electron-phonon interaction of the order 
of X p = 1.3 - 2.3 |]§. 
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FIG. 3. T c versus g for BCS and Nam's model respec- 
tively. 

The BCS parameter A(0)/T c can be easily calculated 
from Eq. © as @ 

A(0)/T c = T d /[T C sinh(l/g)}. (5) 
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FIG. 4. The BCS parameter A(0)/T c versus g. 

As is shown in Fig. [|, A(0)/T c is a function of g or 
T c /Td, and increases with increasing g or T c /Td- In the 
range of g < 0.2 or T c /T d < 0.0076, it has a constant BCS 
value. In fact, this range corresponds to the case of low 
T c superconductors (LTS). In the range of 0.5 < g < 1.5 
it increases almost in a linear of g 7 and has a saturated 
value of 2 for large g. 

In the sprit of Bardeen P0|] , the normal fluid density 
p n (T) = p — p s (T), within the pairing theory, may be 
written as 
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p n (T)/p = 2 / d(w/T)n(w)f(u>/T) [1 - /(w/T)] , (6) 



where f{x) is the usual Fermi function l/[l + exp(x)] and 
the density of states n(uj) — N(ui)/N(0) is given by (|] 

n(u)) = q(uj/T d ) + riBcs(u)r(uj/Td), (7) 
q(w/T d ) = (2/tt) tan-^w/Td), (8) 
r{uj/T d ) = (2/Tr)tan- 1 [n BCS (uj)T d /uj}, (9) 

n BCS (u) = Re{cj/(uj 2 -A 2 ) 1 / 2 }. (10) 

Physically, p n {T) would be resulted from the single parti- 
cle excitation not pairs. Thus, the factor f(x) in Eq. (|^) 
is the occupation probability of the state \k T> and the 
factor [1 — f(x)] is the unoccupation probability of the 
partner state, say, — k i> , and vice versa, respectively. 
The factor 2 comes from the spin sum. The states of 
Eq. (^) are reflections of states being not participated in 
pairings. A word of caution is in order. The u) is the 
dynamical energy which has the kinetic as well as in- 
teraction parts. Physically, the sum of spectral weights 
outside T d < |efc|, result in the states of Eq. (||). Thus, 
the low energy states are realized. In other words, carri- 
ers which do not participate in pairings yield the linear 
T dependence of p n (T) at low T. In fact, these states 
results in the linear T dependence of A(T) at low T. To 
see this, by inserting Eq. (||) into Eq. (||), one can get the 
variation of A(T) at low T as §,|§ 

AA/A(0) = \p n {T)/p= (T/T C )(T C /T d ) (2/tt) In 2, (11) 

similar to the result by d-wave picture ]?]] , 

[AA/A(0)] d = (T/T c ) (T c /A ) In 2 (12) 

via rid (<*0 = <■<->/ Aq, where Ao is the maximum value (anti- 
node) of the order parameter. 

For the quantitative calculations of A(T), we have de- 
termined T c /T d or g [Eq. (|)] via Eq. @, by taking 
the slope of [A(0)/A(T)] 2 near zero temperature. Once g 
or T c /Td is set, no adjustable parameter is used in our 
calculations of Eq. (||) . 

As is shown in Fig. [|, we have obtained good agree- 
ments between calculations and data of BSCCO by Lee 
et al |[ , HBCCO by Panagopoulos et al @ and LSCO 
by Panagopoulos et al ||, and Sr214 by Bonalde et al ||, 
respectively. We picked up not all of data points in the 
papers for clarity. Bonalde et al [|| reported that their 
data at low T vary as T 2 which are resulted from scatter- 
ings by impurities or defects. In a finite Td picture, the 
impurity scatterings make some states at the Fermi level 
not participate in pairings, and result in the T 2 term in 
A(T) at lowT [gT|| . 
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FIG. 5. The temperature dependence of [A(0)/A(T)] 2 
(solid lines) compared with the experimental data for HBCCO 
§, BSCCO ||, LSCO § and Sr214 §. 

However, as shown in Fig. ^, we have obtained poor 
agreement near T c between calculation and data of 
YBCO by Hardy et al |l[] and anisotropic data of YBCO 
by Kamal et al §]. The YBCO b case is good. 
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FIG. 6. [A(0)/A(T)] 2 (solid line) compared with data for 
bulk YBCO [[ij, and a- and b-axes, respectively ||. 

The notion of a finite Td results in the incomplete con- 
densation at zero temperature. By considering the sum 
rule, we can calculate the fraction of states, R, being not 
participated in pairings as [0] 



R{z) 



[nBCs(v) - n{w)} dui/A 



(13) 



A 
A 



n{uj)dw I A 
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(2/n)tan' 1 (z) - (l/zn)ln{l + z 2 ), 



where z = A/T d = [A(Q)/T C ](T C /T d ) which is a function 
of g or T c /T d . 
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FIG. 7. The fraction of states, R, being not participated 
in pairings versus g. 

In Fig. [7] is shown R{z) as a function of g. In the range 
of g < 0.2 or T c /T d < 0.0076, R(z) is negligible. As 
stated before, this range corresponds to the case of LTS. 
Thus, the linear T dependence of A(T) at low T is hardly 
observed in LTS. 

In summary, even though the model of Eq. (||) is ideal, 
the quantitative calculations account very well for data 
for all T ranges without any adjustable parameter, ex- 
cept for YBCO data near T c . Perhaps, the Fermi surface 
effect would play an important role for A(T) in the case of 
YBCO. Of course, the retardation and non-local effects 
should be taken into account as well for improvement. In 
all, the calculations are quite satisfactory and theoreti- 
cally sound. We suggest the pairing interaction energy 
range T d in HTS may be of the order of 1 ~ 2 times T c . 
The linear T dependence of A(T) at low T does not imply 
nodes in the order parameter, contrary to general belief. 

In the spirit of a finite T d , it is recently shown |2l| that 
the spinless impurity scatterings suppress T c and destroy 
superconductivity. Some states at the Fermi level are 
shown not to participate in pairings when there are scat- 
tering centers such as impurities, and result in the linear 
T term in the specific heat at low T. The quantitative 
calculations |E| account well for the reduction of Tc ]23j ] 
and the specific heat data |24|,|2{| in the Zn-doped YBCO, 
respectively. 
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